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1 Introduction 



Dynamical theories exhibiting gauge freedom are described by singular Lagrangians or in the Hamilto- 
nian formalism introduced by Dirac [1,2]. The existence of constraints reduces the true degrees of freedom in 
the system. This reduction has specific features depending upon whether we are in velocity space (Lagrangian 
formalism) or in phase space (Hamiltonian formalism). In particular the dimensions of the constraint sur- 
faces are different; One can prove [3] that the number of constraints in the Hamiltonian formalism equals the 
number of constraints in the Lagrangian formalism plus the number of independent gauge transformations. 
In other words, the number of degrees of freedom in phase space seems to be less than the number in velocity 
space, but getting rid of the gauge freedom in the phase and velocity spaces eliminate spurious degrees of 
freedom. We will show that therefore the number of degrees of freedom are the same in the two formalisms. 

There are two different stages in the reduction of the degrees of freedom. The first is provided by the 
constraints that naturally arise in the formalism by pure consistency requirements. The second corresponds 
to the gauge fixing procedure. In our case it will consist in the introduction of new constraints in such a way 
that the gauge degrees of freedom are thoroughly eliminated. 

In [4] an analysis was carried out on how to implement a gauge fixing procedure in the Hamiltonian 
formalism to obtain the true number of physical degrees of freedom. In the present paper we will proceed 
to general Lagrangian gauge fixing in a way completely independent from the Hamiltonian method. 

We will emphasize the double role played by the gauge fixing procedure. Some constraints are needed 
to fix the time evolution of the gauge system, which was undetermined to a certain extent. The rest of 
the constraints eliminate the spurious degrees of freedom that are still present in the setting of the initial 
conditions of the system. This double role may be relevant for quantization (we treat only classical systems 
here). It has not been adequately emphasized in the literature. 

The gauge fixing procedures in the Hamiltonian and Lagrangian formalisms are technically different. 
For instance, there is no Poisson Bracket available in velocity space for singular Lagrangians, and the 
relationship between velocity and position in the tangent bundle does not carry over to the cotangent bundle. 
This difference is why there are difficulties in implementing a Lagrangian gauge fixing procedure which is 
independent of a pullback of the Hamiltonian one. We here show how to fix these diflaculties and how to 
perform in full generality gauge fixing in both formalisms. Hence our results will include a general proof of 
the matching of the degrees of freedom in Hamiltonian and Lagrangian formalisms for gauge theories. 

Our paper is organized as follows: In section 2 we develop a detailed version of the Hamiltonian gauge 
fixing procedure (with some improvements to [4]). In section 3 we study the Lagrangian gauge fixing proce- 
dure and establish the theorem that the number of degrees of freedom in the Hamiltonian and Lagrangian 
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formalisms are equal. In section 4 we describe some details concerning reparameterization invariant theories, 
since in section 5 we apply our results to two such theories. These examples are the parameter-independent 
Lagrangian for geodesies in special relativity and the Types I and IX spatially homogeneous cosmologies, 
and we will emphasize the role of time-dependent gauge fixing conditions. 

Throughout the paper wc will assume that some regularity conditions arc fulfilled: The Hessian matrix 
of the Lagrangian with respect to the velocities has constant rank, ineffective constraints (such that their 
gradient vanishes on the constraint surface) do not appear, and also the rank of the Poisson Bracket matrix 
of constraints remains constant in the stabilization algorithm (so that a second class constraint can never 
become first class by adding new constraints to the theory). 

We emphasize that we always maintain the equivalence between the Lagrangian and the Hamiltonian 
formalisms [5]. This equivalence holds even before the implementation of the gauge fixing conditions; in 
particular we do not modify the Hamiltonian formalism by adding ad hoc constraint terms as Dirac has 
proposed [2,6,7]. This proposal has been proved to be unnecessary under our regularity conditions [4]. We 
will discuss this point more fully in the Conclusion. 

2 Hamiltonian Gauge Fixing Procedure 

We start with a canonical formalism using Dirac's method, starting from a singular Lagrangian L{q^, q*) 
(i = 1, • • • , N) (q* = dq^ /dt). The functions Pi(q, q) = dL/dq^ arc used to define the Hessian Wij = dpi/dqj, 
a matrix with rank N — P (wc assume this rank is constant), P being the number of primary constraints. The 
Legcndre map from velocity space (tangent bundle for configuration space) TQ to phase space (cotangent 
bundle) T*Q defined by Pi = Pi{q, q) defines a constraint surface of dimension 2N — P. 

The function El ■= Piq^ — i in velocity space (the so-called energy function) is mapped to a function 
on the constraint surface, and in phase space a canonical Hamiltonian He may be defined which agrees with 
this function on the surface. He is not unique, and to it may be added a linear combination X^tpp{q,p) of 
constraint functions •il>p{p= 1, • • • , P), the being arbitrary functions of time t and the vanishing of '4'p{q,p) 
defining the primary constraint surface. These primary constraints may be chosen so that some are first 
class (their Poisson Brackets with all the constraints weakly vanish, that is, vanish on the constraint surface) 
and some are second class (the matrix of their Poisson Brackets with each other is nonsingular — there must 
be an even number of second class constraints or none). 

Let tpl,i>p denote first and second class primary constraints, 2 being the respective A functions. The 
time derivatives of the second class constraints yield 

iPl = {i;l,Hc} + XUi^l,i^l} , 
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equations which may be solved for the A2 by requiring tp^ = 0. These functions A2 are then inserted into 
the expression He + AjVp + •^2'^^ to yield a new candidate for the Hamiltonian: 

He + A?Vp + A^Vp = Hh + AjVp ■ 

The time derivatives of the primary first class constraints yield 

i,l = {i,l,Hc} , 

which can either be zero on the constraint surface or not. In the latter case, new constraints are found. The 
time derivatives of these new constraints will involve their Poisson Brackets with H^ and with the tpj, (and 
will also involve their partial time derivatives). The requirement that these time derivatives vanish will be 
equations, some of which can be solved for some of the Aj functions. The rest of these equations may yield 
more constraints, and the process of requiring the vanishing of their time derivatives is repeated at a deeper 
level. 

This brief description of the stabilization algorithm is not meant to be rigorous, but the process does 
eventually finish. Once the stabilization algorithm has been performed, we end up with [2,5,8]: 

1: A certain number, M , of constraints. These constraints may be more numerous than the ones introduced 
above (M > P), but they arc arranged into first and second classes. The first class constraints have 
weakly vanishing Poisson Brackets with all the constraints, and the matrix of the second class constraint 
Poisson Brackets is nonsingular. These constraints restrict the dynamics to a constraint surface within 
T*Q of dimension 2N - M. 

2: A dynamics (with some gauge arbitrariness) on the constaint surface which is generated, through Poisson 
Brackets, by the so called Dirac Hamiltonian: 

Hn := Hfc + A^.^^, . 

Hpc is the first class Hamiltonian, obtained by adding to the canonical Hamiltonian He pieces linear in 
the primary second class constraints, t/i^ (/U = 1, • • • , Pi) are the primary (hence the superscript 1) first 
class constraints. The secondary and higher first class constraints, obtained from the time derivatives 
of the (f)j^, are not used here. The A** are arbitrary functions of time (or spacetime in field theories). 

3: A certain number (Pi) of independent gauge transformations generated, through Poisson Brackets, 

by functions (/i = 1, • • • , Pi) which have the following form [4,9,10,11]: 

G, = e.^" + e«<A^-^ + e(f + • • • + e^^^-'Ul , 
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where is an arbitrary infinitesimal function of time; e\[' is its r-th time derivative; K/j, is the length 
of the stabilization algorithm for the primary first class constraint and . . . , (j)^" are secondary 
through Ki^-avy, first class constraints. It turns out [12] that one can take these gauge generators in 
such a way that all the first class constraints are involved once and only once in the G^, and so their 
total number equals 

Now we arc ready for the gauge fixing procedure. Even though the order of introducing the gauge fixing 
constraints is irrelevant, we will proceed in the way that makes the whole procedure more illuminating from 
the theoretical point of view. As we said in the Introduction, we will distinguish two different steps in the 
gauge fixing procedure, corresponding to evolutionary laws and initial conditions [4]. In the first step we fix 
the laws of evolution, which otherwise have a certain amount of mathematical arbitrariness. In the second 
step we eliminate the redundancy of initial conditions that are physically equivalent. 

The arbitrariness in the dynamics is represented by the Pi functions A**. To get rid of this arbitrariness, 
we introduce a set of Pi constraints xji — (m = li " " " ) -Pi), defined so that their own stability equations, 
under dynamical evolution, will determine the functions A**. To this end we must require that the matrix 

■■= {xl<l>l} 

be non-singular. The conservation in time of this new set of constraints leads to 

= ^ = ^ + ^^^^ + "^^^ ' 

which determines as 

X'^ = _^C-T{{xl,Hpc} + ^) . 

The dynamical evolution thus becomes completely determined. The imposition of these constraints causes 
the dynamics to be further restricted to the (2A/' — M — Pi)-dimensional constraint surface defined by xji = 0- 

The gauge fixing procedure is not yet finished. It is necessary to address the issue of initial conditions, 
which we call "point gauge equivalence," our second step. Let us clarify this crucial point: Even though the 
dynamics has now been fixed, there is still the possibility of gauge transformations which take one trajectory 
into another. To check whether these gauge transformations do exist, we need only check their action at a 
specified time. That is, the points on the set of trajectories at a specified time are unique initial data for 
the trajectories. If a gauge transformation exists which relates two initial data points, then these two points 
are physically equivalent. We will obtain the generators of the transformations which take initial points into 
equivalent ones ("point gauge transformations") and use them to fix the gauge finally. 
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Consider the gauge generators given above at, say for simplicity, t = 0: Gi_i{0). The most arbitrary point 
gauge transformation at f = will be generated by G(0) = X^^L^ ^/i(0)- The arbitrary functions and 
their derivatives become, at the given time t = 0, (infinitesimal) independent arbitrary parameters (there 
are F in number). We redefine them as 

elf ^"'"^0) ; (/z = 1, • • • , Pi) (z^ = 1, • • • , i^^) . 



These point gauge transformation generators must be consistent with the new constraints xji- This 
requirement introduces relations among the a^^,^ : 

Pi K„ Pi 

Remember that the matrix C^^ = {xlt't'ji} is nonsingular. These relations imply 

Pi 

ap,i = -(C-i)'"^^^a^,.Jxi,</';^}- 
As a consequence, the independent point gauge generators are 

4>ir ■■= <i>ir - <i>lic-'r{xl, (m = 1, ••• , Pi) (v = 2, • • • , K,) . 

Notice that new point gauge generators only exist when there are secondary first-class constraints, that is, 
when the length of at least one of the is greater than one. 

Recall that F is the number of first class constraints in the original theory, including primary, secondary, 
and higher constraints; we conclude that there are F — Pi generators that relate physically equivalent initial 
conditions. To eliminate the extraneous variables, we will select a unique representative of each equivalence 
class by introducing a new set of F — Pi gauge fixing constraints, x^" ~ (/i = 1, • • • , Pi) (i^ = 2, • • • , Kfj,), 
such that 

det|{x;-,e}|7^0, v^l, j^j^l, 
in order to prevent any motion generated by (pl" . The stability requirement is 

which ~ means vanishing on the constraint surface; this requirement simply dictates how the x^" evolve 
off the initial data surface. 

Notice that we have explicitly allowed time dependence in the x^" constraints. In fact, time dependence 
is necessary in the special case when Hpc is a constraint (first class, of course). This point will be clarified 
in Section 4. 



This ends the gauge fixing procedure. Now we can count the physical number of degrees of freedom: 
The M constraints left after the stabilization algorithm restricted motion to a 2A'' — M-dimensional surface 
in T*Q. The gauge fixing constraints needed to fix the evolutionary equations number Pi. Finally there 
are F — Pi point gauge fixing constraints needed to select physically inequivalent initial points. (The total 
number of gauge fixing constraints equals the number F of first class constraints in the original theory.) The 
final number of degrees of freedom is 2A'' — M — F. Notice that M — F is the original number of second class 
constraints and is therefore even. Consequently, 2A'' — M — F is even; this result agrees with the fact that 
the above procedure makes all constraints into second class ones, and in this case the constraint surface is 
symplectic [13]. 

3 Lagrangian Gauge Fixing Procedure 

We first use a stabilization algorithm similar to the one used in the Hamiltonian formalism. (In the equa- 
tions below, we use the summation convention for configuration space indices i = 1, • • • , N.) The equations 
of motion obtained from the Lagrangian L are (assuming for simplicity no explicit time dependence): 

where 

'■^ dq^dq^ ' ' dcfdq^ dq^ 
If Wij is singular, it possesses P null vectors 7^, giving up to P (these relations may not be independent) 

constraints 

ai7p - . 

It is easily shown that there exists at least one M*^ and 7^'' such that 
and therefore [14] 

where fj'' are arbitrary functions of t. 



if = M^'as + fjPYp (with 77^ = 7f 9*) 



The stabilization algorithm starts by demanding that time evolution preserve the aj7* constraints. 
Sometimes new constraints are found; sometimes some of the f]'' are determined; eventually the dynamics is 
described by a vector field in velocity space 

the a* are determined from the equations of motion and the stabilization algorithm; rj'^ {iJ- = 1, • ' ' ) -Pi) are 
arbitrary functions of time; and 



where ^^^7^ are a subset of the null vectors of W^, corresponding to the primary first class constraints found 
in the Hamiltonian formalism. It is not necessary to use the Hamiltonian technique to find the F^, but it 
does facilitate the calculation: 

dpi ' 

where the (f)j^ are the primary first class constraints, and Pi stands for the Lagrangian definition of the 
momenta Pi = dL/dq^. The Pi number of t]'' is the same number as in the Hamiltonian formalism. There 
are left M — Pi constraints [3]. 

At this point it is useful to appeal to the Hamiltonian formalism for the computation of the Pi inde- 
pendent gauge transformations. The result includes the definitions of the functions ^jH' , and then in the 
Lagrangian formalism we define 

These functions give the infinitesimal Lagrangian gauge transformations as 

"til V V-Jm ' 

the C/j, being arbitrary functions of time. 

As wc did in the Hamiltonian formalism, the first step in the gauge fixing procedure will be to fix the 
dynamics to determine the arbitrary functions 77^. To this end we introduce Pi constraints, ~ 0, such 
tliat _D^^ T^fii^l has non-zero determinant: det|r^cijj^| ~ det |-D^i^| ^ 0. Then the functions become 
determined by requiring the stability of these new constraints: 

= Xa;° = Xoa;°+r?''r^a;° . 

This relation gives 

which determines the dynamics as 

XF = Xo-{Xou,l){D-^y''r^ . 

Although the time evolution is fixed, as in the previous section there still remain some point gauge 
transformations in the constraint surface that we should get rid of. Again, these transformations may be 
thought as affecting the space of initial conditions. In fact, we can extract those transformations at t = 
that preserve the gauge fixing constraints a;° ~ from the general gauge transformations. This general 
transformation is 
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where now 



dt ^ ^ dt ' 

In this expression we use the values at t = 0: 



5* 



after redefining a^^i^ :— e^^" '"'(0) (we have defined i^-^+i = 0). 



Notice that now, due to the presence of the time derivative of ^^g', runs from to i^^. This is a key 
difference with respect to the Hamiltonian case, where the a^,,^ parameters had indices running from 1 
to Kfjt. We call the result for the independent point gauge transformations (at t = Q) ^(o)wj^: 

/ /r)i fi Pi, ,0 Pi, fi \ \ 

= E E (tf /^v + W^-fU) + if + -.,or,A > 

where ^ = 7^. At this point, recalling that det |r^(i;°| ^ 0, wc sec that the stability conditions ^(o)'^° = 
allow the determination of q^,o in terms of Q:^,i^ (m = 1, • • • , -Pi) (v = !> • ' ' > -^^/i)- 

We conclude that the independent point gauge transformations (5(o) that still remain, relating physically 
equivalent initial conditions, are parameterized by Q:/i,i^ (a* = 1) ' ' ' : -Pi) (v = 1) • • • > -?^^i)- Their number 
equals F, the total number of first class constraints in the Hamiltonian theory. To eliminate these trans- 
formations we introduce F new gauge fixing constraints loIT — (/x = 1, • • • , Pi) (i^ = 1, • • • , K^), with the 
conditions: 

1: The system 5(o)t^^'' = 0, which is linear in the {fJ- = l,---,Pi) (i^ = 1,---,K'^) has only the 

solution a^^i^ = (so that no point gauge transformations are left). 

2: Xp{wlt) ~ (the requirement of stability under evolution). 

Now we have completed the gauge fixing procedure. For reasons similar to the ones raised in the Hamil- 
tonian formalism, there are cases where a time dependent constraint shows up necessarily. Our examples in 
section 5 will be two of these cases, and we discuss these cases in the next section. 

Summing up, the gauge fixing constraints introduced in velocity space (that is, in the Lagrangian 
formalism) are w^" {/i = 1, - ■ ■ , Pi) (i^ = 0, • • • , K^). Their number is P + Pi, and therefore the total number 
of constraints becomes (-M — Pi) + (P + Pi) = M+F. The number of degrees of freedom is then 2N — M — F. 
Comparison with the results of the previous section shows that we have proved 
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Theorem. The number of physical degrees of freedom in constrained Hamiltonian and Lagrangian for- 
malisms is the same. 



Observe that this result, which was obviously expected on physical grounds, is nontrivial. In fact, before 
introducing the gauge fixing constraints, the dimensions of the constraint surface were different in the two 
formalisms. This means that the gauge fixing procedure has to make up for this difference — and we see that 
it does. 

4 Reparameterization Invariant Theories 

Reparameterization invariant theories provide interesting cases for the application of the preceeding 
sections. Examples of this kind, including spatially homogeneous cosmologies of Types I and IX, will be 
treated in the next section. 

If we consider the infinitesimal reparameterization t ^ t' = t — e{t), with e an arbitrary (infinitesimal) 
function, the trajectories q^{t) (any trajectory, not necessarily solutions of the equations of motion) change 
accordingly, q^{t) — > q'^{t'). If we define the functional infinitesimal transformation 5q^ = q'^it) — q^{t), the 
transformations we will consider are of the type 5q^ = eq^ + /i% where the /i' terms involve time derivatives 
of e. The theory is reparameterization-invariant when the Lagrangian remains form-invariant under these 
changes: 

L{q{t),q{t))dt = L{q'{t'),q'{t'))dt' . 

In such a case, we define 

dL ^ j dL ^ ■ 

•= TTl^l' + jTT-Sq' . 
oq^ aq^ 

It is easy to check that SL is a total derivative: 
This equality can be transformed into 

[LW + = , 

where [L]i are the Euler-Lagrange derivatives and where G is the conserved quantity of Noether's theorem: 

G can be expanded as a sum of an e term and terms involving derivatives of e. Each of these terms must 
be a constraint of the theory because G is a constant of motion for whatever arbitrary values we give to e 
(this is a general argument for gauge symmetries). Using the form of of 5q^ introduced above, we obtain, for 
the e coefficient of G, 
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namely the Lagrangian energy function El. Its corresponding canonical quantity is the canonical Hamilto- 
nian He; therefore we have stated the following: 

Theorem. The canonical Hamiltonian (if it is non-zero) in a reparameterization-invariant theory is a con- 
straint. 

In general He will be a secondary, first class constraint, but a particular case is worth mentioning: 
When all the configuration space variables transform as 6q^ = eif ("scalars"), then G becomes G = cEl- 
The relation [L]iSq' + = gives 

[L]ieq' + Cj^El + cEl = . 

Since this relation is identically zero for any function e, we conclude that the coefficient of e, namely El, is 
identically zero. Therefore the Lagrangian is homogeneous of first degree in the velocities: The canonical 
Hamiltonian vanishes in this case. The remaining pieces tell us that is a null vector of the Hessian matrix 
of the Lagrangian, and in case this is the only null vector, that there are no Lagrangian constraints. This 
situation occurs exactly in the case of the relativistic free particle, which is described by the Lagrangian 
L = ^yxi^Xfj, and which will be treated in the next section. 

Now we will prove another result for reparameterization-invariant theories, the need for time dependence 
in some gauge fixing constraint. Suppose that He vanishes. Then to fix the dynamics (that is, to determine 
functions in Hu which do not all vanish) by using the conditions = 0, it is necessary that at least 
one of the constraints have explicit time dependence. If Hq doesn't vanish, then it is a constraint; the 
first class Hamiltonian, Hpc, will necessarily be a first class constraint of the original theory. After the first 
step of the gauge fixing procedure (in which the dynamics is determined) , the final Hamiltonian Ho will be 
a first class constraint that generates motions tangent to the first-step gauge fixing surface. 

This latter result means that Hd will become a part of G(0), the generator of point gauge transfor- 
mations which relate physically equivalent initial conditions. Then, in order to fulfill the two requirements 
introduced in the second step of the gauge fixing procedure, it is mandatory that at least one of the gauge 
fixing constraints be time-dependent: Otherwise there is no way to satisfy the gauge fixing conditions. By 
choosing variables appropriately, we can always end up with only one time-dependent gauge fixing constraint. 
Therefore we have proved: 

Theorem. Reparameterization-invariant theories necessarily require that one of the gauge fixing constraints 
be time-dependent. 

This result is clearly expected from the physical interpretation of this kind of theory: The existence of 
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reparameterization invariance as a gauge symmetry implies that the evolution parameter — the "time" — is 
an unphysical variable. 

5 Examples 

We will first discuss the case of a relativistic free particle and then spatially homogeneous cosmological 
models of Bianchi Types I and IX. 

One Lagrangian for a free particle in special relativity is (r is the path parameter) 

where ' means d/dr. (For convenience we take rj^i, = diag(l, — 1, — 1, — 1).) The action integral J Ldr is 
invariant under arbitrary reparameterizations, so r is not necessarily proper time. The conjugate momenta 
functions in velocity space are 

Tjni/X 

and therefore the velocity-space energy function is 

El = Pm*" -L = 0. 

First, we examine this system from the Lagrangian point of view. The definition of implies 
The equations of motion imply that =const. The dynamics vector is 

To fix the dynamics (to determine the arbitrary function A), we use the constraint that the path parameter 
is proper time: 

this constraint implies A = 0. 

We must now set initial data, which will be seen to be equivalent to fixing the zero point of proper time. 
Because the Lagrangian is homogeneous of first degree in the velocities, the gauge transformation is of the 
form 5x^ = ex^. At r = 0, this becomes 

5ox^ = ax^ , Sox^" = ^3x^' , 
12 



where a, (3 are infinitesimal constants. To be compatible with the proper time constraint, we must have 
(3 = 0. We will choose a constraint % such that Sqx = implies a = (to prevent any point gauge 
transformation); therefore 

The evolution of x obeys 

dx>^ dr ' 

and so x must be explicitly time-dependent. One convenient choice is 

X = — x^T . 

The zero point of the path parameter is set by this requirement: a;'^(0) = 0. There are then a six-parameter 
set of paths which are solutions, the six parameters being the three positions and three spatial components 
of the velocity at r = 0. 

We now turn to the Hamiltonian discussion. The Legendre map to phase space is a map onto the 
surface defined by ij'^'^Pf^pi, = 1. The canonical Hamiltonian He is a function on the surface, but a trivial 
one: He = 0. The actual Hamiltonian in the Dirac procedure is the addition of an appropriate function of 
the constraint to He, namely 

Hd = ^xiTM^P^.Pu - 1) . 

The equations of motion in phase space are p^ = and x^^ = Xiqi^^Pv The one primary constraint (M = 1) 

<k-=\{yi'"'p^.p.-i) 

is first class, and there are no secondary constraints (Pi = 1). The gauge transformations are generated by 
the one function 

G = e{t)(l> . 

The gauge fixing procedure has two steps. The arbitrariness in the equations of motion represented by 
A requires a constraint function x defined so that 

Clearly this function must be time dependent, or else the result of requiring x = will be A = 0. One choice 
is 

X = r-arO ^ 

Po 

The Hamiltonian is now definite: 

H=^{V'"'P^.P.-1) . 
^Po 
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The second step is to consider G(0) = acf), where a = e(0), as a generator of point gauge transformations. In 
order for the Poisson Bracket of x with G(0) to vanish, clearly a must be zero, so that there are no further 
gauge fixing steps to perform. 

The result of gauge fixing in the Hamiltonian treatment of this free particle is therefore that motion in 
phase space is within the surface defined by 'q^'^Pij.Pv = 1 (a 7-dimensional space), and the motion is given 

by 

„ ri^^Pv 
Pa = const , = T + const . 

Po 

Note that there are six free parameters for these paths, since the x constraint fixes the motion to be 
x^ = T. Of course, this is the same number as found in the Lagrangian treatment, though here it was 
convenient to normalize r by the requirement x^ = 1, rather than by the requirement that r be proper time 
(to which it is proportional, anyway). 

Our next set of examples are spacetimc metrics which are invariant under a three-dimensional isomctry 
group which is transitive on spacclikc three surfaces [15]. The metric is best expressed in a basis of differential 
forms which is invariant under the group. One such basis consists of the four one-forms {dt,LO^}, where 

where the C]^ are the structure constants of the Lie algebra of the group and obey 

The second relation is the Jacobi identity, and in the case of a three-dimensional Lie algebra is exactly 
equivalent to 

In what follows, Greek indices range over 0,1,2,3, while Latin indices range over 1,2,3. The summation 
convention will be followed even if both indices are superscripts or subscripts. 

In this basis, the line element is 

ds^ = -N^dt^ + gst{N'dt + uj'){N*dt + w*) , 

where A'' is the lapse function and is the shift vector; N, N^, and Qst are functions only of the time t. 
The models are classified according to a standard listing of the possible structure coefficients into nine basic 
classes, called Bianchi Types. 

Our examples will be the models of Type I (Cjj, = 0) and Type IX (Cjj, = kcijk, Cijk is the completely 
antisymmetric Levi-Civita symbol defined by ei23 = 1). The constant k in the latter models is redundant but 
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is included to allow the limit A; = 0. In both cases the shift vector may be taken to be zero and the spatial 
metric Qst may be taken to be diagonal (we treat only vacuum models). In both cases (as well as in any 
model which has = 0), the Lagrangian of the system may be calculated from the spatially homogeneous 
form of the scalar curvature R. 

Instead of the above basis, it is more convenient to use an orthonormal basis {<J^}- Since we take = 
and gij diagonal, this basis is defined by 

cr° = Ndt , cr* = e~^e'^*a;' (no sum on i) , 

where O and Pi are functions only of t, with J^Pi =0. The metric components Qij are given by 

(ff,,) = cliag(e-2^+2/'0- 

In this basis the line element is 

ds"^ = ■q^.^a'^a'' , r?^^ = diag(-l, 1, 1, 1) . 

Since J20i = 0' define /3± by 

Pi=P+ + \/3/3_ , p2=P+- \/3/3_ , ps = -2/3+ . 

It is also possible to reparameterize time to make A?' = 1 (or some other function), but it will be seen 
that N must remain a dynamical variable, at least at first, in order not to spoil the Lagrangian procedure. 

The field equations for the functions N,fl,f3i are the Einstein equations (for a vacuum). The Ricci tensor 
coefficients are functions only of t in this case, so the equations are ordinary differential equations. The 
most convenient form of these equations for our purposes will involve the Einstein tensor G^i, = R^u — ^Rilnv 
(The equations may be written either as G^;^ = or R^^ = 0.) 

Our first example is Bianchi Type I, in which Cjj, ~ 0. In this case the one-forms are expressible in 
terms of coordinates as = dx^. The line element is 

ds^ = ri^^a'^a" , a° = Ndt , a* = e'^e'^'dx^ (no sum on i) . 

The appropriate Einstein tensor components (or rather, independent linear combinations of them), which 
are to be set equal to zero, are 

^(Gn + G22 - 2G33) = /3+ - ^/3+ - 3/3+0 = , 

iV2 .. N . . . 

= (Gn - G22) = /?- - ^/3- - 3/3_0 = , 



^(Gn + G22 + G33 + 3Goo) = O - - 30^ = . 
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The action integral for general relativity is X = J R^/\g\d'^x (up to an irrelevant proportionality con- 
stant), where R is the Ricci scalar and g is the determinant of the metric in coordinates {x^}. In our case 
we have 

J = j Ra° Aa^ Aa^ Aa^ = j RNe'^^d^x . 

The result here is 

We take the volume of space / d^x = 1/12 and integrate by parts with respect to time, dropping the endpoint 
contributions: 

Variations of I with respect to N, (3+, give the Einstein equations listed above. 

Note that reparameterization of time can be used to set A'' = 1 (or some other function), but if N is 
eliminated from the Lagrangian, the Gqo = field equation (which is a constraint equation) will not be 
derivable. 

But of course the Lagrangian L, 
is singular, since N doesn't appear: 

dL „ 
Pn = — - = . 

dN 

An attempt to form the Hamiltonian thus yields 

He is a function in T*Q, namely in principle a fimction of N, and pN,pn,P+,P-, but it happens 

to be independent oi pn- The Legendre map from TQ to T*Q is 

{N,n,p+,p_} = {N,n,p+,p_} , 

PN = ,PQ = ,P+ = -^/^+ ,P- = ~/v~^" ' 

and thus maps the 8-dimcnsional TQ into the 7-dimensional subspace of T*Q defined hy pn = 0. We thus 
identify the primary constraint in T*Q: 

Pat ~ , 

where ~ means that every solution of the equations of motion has to satisfy the constraint. Since there is 
only the one primary constraint, it is first class. 
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Therefore, to He may be added an arbitrary function of pn which vanishes when pj^ = 0, the simplest 
being p^ itself: 

H = Hc + XpN ■ 

This arbitrariness means that we have a gauge- type freedom. The time derivative oi pn is given by 

= {p^,H} = = -he'^'i-ph +pI +p'-) . 

The requirement that pjv — thus implies that Hc/N ~ or —p^ -\- p\ + p^_ — Q. 

Now is the time to generalize this example, so as not to repeat a lot of material. The generalization will 
still be quite concrete: The Type IX cosmology has 

Cjfc = keijk ■ 

In this case, the action integral is 



J = j Ra° Aa'^ Aa^ Aa^ 



where the orthonormal frame of one-forms is 

cr° = Ndt ; cr* = e~'^e'^'D,'- (no sum on i) ■ dw"- = \keijk^^ A uj'^ . 

Thus I is 

J = y" RNe-^^ dtAw^AiJ-Auj^, 
and the spatial integral may be set equal to 1/12 as before. 

The result for the Lagrangian is are defined as before): 

-an 

i = ^(-fi' + /3^ + /3i)-fc'Are-M/3+,/3_) , 

where the function V is 

= \e-^^+ - 2e-2^+ cosh(2V3/3_) + e'*''+ (cosh(4 V3/3_ - 1) . 
The particular form of V is not important, nor are the specific forms of the field equations. 
Needless to say L is still singular: 

Pjv = . 

There is still just one primary, first class constraint, p^r ~ 0, and the Hamiltonian is given by 

H = Hc + XpN , 
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where 

He = iNe'^i-pl +p\ +p2 ) + fc2Are-^l/(/3+,/3_) . 
Stability of the primary constraint gives 

Pn = {pn, He} = — ^ ~ , 

where tp is the secondary constraint 

V = h'^'ipl + pi - pI) + fc'e-^y (/3+, . 

At this point we can see that the canonical Hamiltonian is just a constraint: He = Nil). No further 
constraints appear. 

Primary and secondary constraints are both first class, the gauge generator is made out of them: 

G = epN + ■ 

After the redefinition of the arbitrary function, e = Nr], the transformations are 

(5/3+ = /?+?7 , 513- = $-Tj ,6n = nr] ,SN = Nr] + Nri . 

(Thus, /3+,/3_,f2 behave as scalars under reparameterizations, and A'' behaves as a vector.) 

According to the theory developed in the previous sections, there will be two gauge fixing constraints, 
one of them time-dependent. These two constraints have to make the two original constraints second class 
in the theory in order to determine the arbitrary function in the DLrac Hamiltonian. The simplest way to 
proceed is to reverse our formal methodology and first to write down the initial data gauge fixing constraint. 
If this constraint does not depend on N, its stability will give us a new constraint which is AT-dependent. 
Then the stabilization of this new constraint will fix the dynamics. 

We start with the time-dependent constraint (a simple and common choice): 

X^^^ :=f2-t~0 . 

Its stability leads to the requirement 

X(i) = {n, Hc}-1 = -iVe-3 V - 1 ~ . 
Prom this expression we get the new gauge fixing constraint 

.= TV + — ~ . 

Pn 
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Notice that since N is required to be positive (on physical grounds), this constraint imphes pa <0. 



In its turn, stability of x^^\ because of the presence of N, will determine the arbitrary function A in 
the Dirac Hamiltonian Hjj = He + XpN- Actually, in order to fix the gauge in the equations of motion, we 
do not need to know the value of A: It only matters for the equations of motion for the variable N, and 
the constraint x^^^ ~ already relates N to other variables. Therefore, since we are concerned with the 
evolution of /3+, /?_, and their associated momenta, the Hamiltonian is simply: 

H = -^(pI +p'- -p^ + 2fc2e-^^y(/3+,/3_)) , 

where we have used the constraint x''2) — to substitute for N. This procedure is correct because N was 
an overall factor of a constraint in He- He = N^p. This result can be rephrased as follows: We have gotten 
rid of a couple of canonical variables, N,pN, by using the Dirac Bracket. In this case, the Dirac Bracket for 
the rest of the variables is the usual Poisson Bracket.. 

The constraint ^ ~ can be used to define the evolution of po, in terms of the other variables (since the 
constraint Q. = t fixes the time dependence of Q). It is conveniently factorized as 

i^ = -\e'''(pn + HR){pn-HR) , 

where 

Hr := ^p\+p'_+2k^e-^^V{f3+,(5-) . 

This factorization implies 

H=^{pa-HR){pn + HR) . 
2pn 

Satisfaction of the constraint -0^0 implies + Hr ~ (the physical interpretation of the model 
requires pn < 0). Then pQ — Hr ~ 2pa, and the Hamiltonian can be equivalently written as 

H = pa + HR. 

Observe that this constraint defines the evolution of pa in terms of the other variables. Thus if we restrict 
ourselves to the evolution of the variables /?_ and their canonical conjugates p+,p-, then the dynamics in 
this reduced space is described by the Hamiltonian Hr. Once the evolution of these variables is determined, 
the rest of the variables (four in number) have a time-evolution dictated by the constraints (also four in 
number) of the theory. 

With regard to the Lagrangian formulation, there is only one constraint in velocity space: It is the 
puUback of the secondary constraint ip, since the puUback of the primary one is identically zero. The 
dynamics (given by a vector field in velocity space) has only one arbitrary function, which multiplies the 
vector field T = d/dN. 
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Let us briefly sketch the Lagrangian gauge fixing procedure in this case. If we start again with the time 
dependent gauge fixing constraint — t ~ 0, its stabiHzation will require O — 1 ~ 0. The stabilization of this 
new gauge fixing constraint will give a new constraint which is dependent on N (plus other velocity space 
variables). Finally, the stabilization of this last gauge fixing constraint determines the arbitrary function of 
the dynamics. We end up with three gauge fixing constraints that add to the original single constraint of the 
theory to give the elimination of four degrees of freedom. This numbering is the same as in the Hamiltonian 
case. 

6 Conclusions 

A singular Lagrangian, that is with singular Hessian Wij = d^L/dq'^dq^ , results in constraints, a dy- 
namics with some arbitrariness, and gauge transformations which reflect this arbitrariness. The Legendre 
map to phase space defined by = dL/d(f therefore maps the 2N dimensional velocity space TQ to a lower 
dimensional surface in phase space T*Q. The requirement that the dynamics on this surface be consistent can 
be used to reduce its dimensionality somewhat, but there still remain the same three ingredients: constraints, 
arbitrariness in the dynamics, and gauge transformations. Consistency requirements in the Lagrangian for- 
malism can also be used to reduce somewhat the dimensionality of the constraint surface in velocity space, 
but in general the constraint surface in phase space will have smaller dimensionality than the constraint 
surface in velocity space. The number of gauge transformations in TQ and in T*Q is the same. 

In this paper we show how to determine the dynamics and fix the gauge in both the Hamiltonian and 
Lagrangian formalisms. We believe the Lagrangian discussion is new and useful. The result includes a proof 
that the final number of degrees of freedom in the two formalisms is the same. Important parts of our 
methods are the two steps of determining the dynamics and determining the independent initial data. 

The first step in both the Lagrangian and Hamiltonian formalisms is the stabilization algorithm (which 
ensures consistent dynamics). The result is M constraints in the Hamiltonian case (a constraint surface of 
dimension 2iV— M in T*Q) and M — Pi constraints in the Lagrangian case (a constraint surface of dimension 
2A'' — M + Pi in TQ); Pi is the number of primary, first class constraints. 

There are Pi gauge transformations (in both the Lagrangian and Hamiltonian cases), and they are 
formed using a total of F functions (first class constraints in the Hamiltonian formalism). In the second 
step, in both the Lagrangian and Hamiltonian cases. Pi gauge-fixing functions are used to determine the 
dynamics (to obtain the dynamics vector Xp or the Hamiltonian Hjj respectively). The dimensionality of 
the constraint surfaces in each case is reduced by Pi . 

The third and final step is to use gauge fixing functions to determine physically inequivalent initial 
data. In the Hamiltonian case F — Pi functions are used, so that the number of degrees of freedom is 
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2A'' — M — Pi — {F — Pi) = 2N — M — F. In the Lagrangian case F functions are used, and again the number 
of degrees of freedom is 2A'' — M + Pi — P\ — F = 2N — M — F. The Hamiltonian method directly shows 
that 2N — M — F is an even number. 

We have emphasized the double role played by the gauge fixing procedure. Some constraints arc needed 
to fix the time evolution of the gauge system, which was undetermined to a certain extent. The rest of 
the constraints eliminate the spurious degrees of freedom that arc still present in the setting of the initial 
conditions of the system. This double role has not been adequately emphasized in the literature. Although 
we only treat classical systems in this paper, our approach should also be relevant to quantum ones. 

Let us briefly comment on Dirac's extended Hamiltonian formulation [2]. Dirac suggests that the 
canonical dynamics be modified by adding all secondary first class constraints to the Hamiltonian in an ad 
hoc manner. The result is as many arbitrary functions in the dynamics as first class constraints, and every 
first class constraint generates an independent gauge transformation. In our dynamics fixing step, we would 
introduce F constraints, and the process would then be finished (the initial data fixing step would be empty). 
The counting of degrees of freedom still agrees with ours. Furthermore, if we take one of our sets of gauge 
fixing constraints x^T — 0, then this set works for the extended formalism also, and we end up with the same 
dynamics. However, one can also use a more general gauge fixing procedure for the extended Hamiltonian 
theory (without the stabilization condition we required for the X^"); in this case the extended Hamiltonian 
theory will not necessarily be equivalent to the Lagrangian formalism. 

We have discussed, also, the special case of reparameterization-independent theories. In those cases, as 
we showed, the canonical Hamiltonian He is a constraint of the motion. In general He will be a secondary, 
first class constraint, but in some cases it will vanish. In reparameterization-independent theories, we showed 
that one of the gauge fixing constraints must be time-dependent. In these theories, the "time" is defined 
through some dynamical variables of the system. Since there is a good deal of choice for this definition 
of time, one has to be very careful in verifying that the time-dependent constraint is consistent with the 
physical interpretation of the model. 

We have illustrated these ideas by considering two interesting cases. The first is the motion of a free 
particle in special relativity using the Lagrangian L = ■^ry^^i^i". This Lagrangian is reparamcterization- 
invariant, and the canonical Hamiltonian vanishes. The result of applying our procedures is that the paths 
arc parameterized by six parameters (three positions in space and three velocity or three momentum initial 
values), and that time may be parameterized by proper time, as one does expect. 

The second case included two spatially homogeneous cosmological models in general relativity, the 
vacuum Bianchi Type I and Type IX models. In these models, the number of degrees of freedom is found to 
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be reducible to four, and the role of the time-dependent gauge fixing procedure is clarified. 

We have previously [16] described some of the above results and intend to extend many of these ideas, 
for instance by looking into aspects of field theory. For example, the various gauges used in electromagnetic 
theory, including the Coulomb (V • ^ = 0), Lorcntz ~ 0), and radiation {Aq = 0) gauges, apply in 

different formulations, cither Lagrangian or Hamiltonian. How they affect the true degrees of freedom of the 
electromagnetic field may be clarified by our methods. 
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